Effect of Multi-Year Environmental and Meteorological Factors on the Quality Traits of Winter Durum Wheat

A detailed study was made of the effect of rainfall, average temperature and hot days on the gluten index and Minolta b* value of winter durum wheat sown in the field in 16 consecutive crop years (2005–2020). The joint analysis of these two technological quality traits represented a complex (plant–environment–meteorological factors) approach for the identification of durum wheat cultivars carrying an optimum combination of the two traits and for the determination of quality stability. The results of GGE-biplot analysis indicated that the cultivar that had the most favorable combination of the traits was ‘MVP’, while cultivar ‘GKS’ had the best gluten strength and ‘MVH’ the best yellow pigment content. Correlation analysis and stepwise regression between various meteorological factors (rainfall, mean temperature, number of heat days per 10-day period during grain-filling) and the two technological quality traits indicated that the expected value of the quality traits could be reliably estimated based on meteorological factors, with a generally negative effect on gluten index and a positive one on yellowness in all cultivars.


Introduction
The vast majority of the durum wheat produced and marketed in the world is grown under non-irrigated, semi-arid conditions [1].
The main growing areas are the North and South Mediterranean zones and the Lusitanian zone [2]. Three methods are used in Europe for the production of durum wheat [3]. In countries in Southern Europe spring cultivars are sown in autumn. Due to the mild winters, the crops are unlikely to suffer cold damage, while they are able to utilize the greater quantity of precipitation falling in winter. In cooler, more northerly areas (Austria, Germany) the cultivars sown are again mostly spring types, but as the risk of freezing is greater in the colder winters, farmers sow the seed in spring. In Central and Eastern Europe (Pannonian region) most of the cultivars grown are true winter durum wheat that is sown in autumn. The main reason for this is that in this region, winter types are able to utilize the better quantities of rainfall available in autumn and winter [4], so their yield potential is substantially better than that of spring-sown crops [5,6].
Winter durum wheat cultivars have also been registered in the United States [7] and Canada [8], but in these countries, the ratio of winter cultivars grown is negligible compared with the spring genotypes.
Durum wheat stands are in part exposed to diverse abiotic stress factors in different environmental zones. While drought and heat stress influence plant development in the North and South Mediterranean regions, in the Pannonian zone these are combined with cold stress. The strong gluten matrix and the high content of yellow pigment are two of the most important components of the technological quality of durum wheat and both present an important breeding target for this crop [9]. While the former plays a role in determining the cooking properties of its main processing product, pasta, the latter is important in heading for a period of 15 days, the bread wheat cultivar investigated had a significantly lower Zeleny sedimentation value, but in the same experiment, there was a significant increase in this trait for the winter durum wheat cultivar 'Mv Makaróni' [29].
Based on this, in the present work, five winter durum wheat cultivars were sown under field conditions between 2005 and 2020 in order to carry out a complex analysis of plant-environment-meteorological factors aimed at obtaining an accurate picture of the magnitude of changes in individual parameters and of the effect of environmental and meteorological factors on the quality traits of winter durum wheat. The analysis of performance using data from multi-year trials allowed a comprehensive, comparative understanding of the germplasm evaluated and of the crop years when the trials were conducted [30]. The specific aims were (1) to explore the effect of the year on the gluten index and Minolta b* value of winter durum wheat cultivars, (2) to evaluate the genotypes for both quality traits and their stability across environments, (3) to identify meteorological factors with a decisive effect on gluten index and Minolta b* value, and (4) to analyze the combined effect of meteorological factors on the gluten index and Minolta b* values of the cultivars tested.

Effect of the Genotype, the Crop Years and Their Interaction on the Gluten Index and Minolta b* Values of Winter Durum Wheat Cultivars, Based on 16 Consecutive Crop Years
In the statistical analysis, significant differences were demonstrated between years and durum wheat genotypes on the basis of the critical values for F distribution at the relevant probability levels and degrees of freedom. ANOVA analysis showed the significant effect of the crop years, the genotype and their interaction on gluten index and Minolta b* values, based on the critical values for F distribution at the relevant probability levels and degrees of freedom (Table 1). The influence of the crop years on the gluten index and Minolta b* values of five winter durum wheat cultivars ('GKB', 'GKS', 'MVM', 'MVP' and 'MVH') grown in 16 consecutive years (2005-2020) was shown in Table 2, leading to the identification of the meteorological factors having a decisive effect on these traits. Judging from Tukey's HSD test, the 2008/2009 growing season was the most favorable for the development of a strong gluten matrix (mean gluten index value = 89.81) ( Table 2). The lowest gluten index was recorded in 2019 (27.48), though the value obtained in 2012 did not differ significantly from this (34.00). Large significant differences were found between the genotypes: "GKS" and "MVP" had the highest and most stable gluten index values, while the "MVM" variety had the lowest in each crop year. The Minolta b* values ranged from 21.42 to 28.14 in the different crop years, the least favorable crop year was clearly seen to be 2017 (21.42). All the cultivars with the exception of 'GKS' (24.18) and 'MVM' (24.39) had significantly different Minolta b* values. 'MVH' was ranked first, with an outstandingly high value (27.19), followed by 'MVP' (24.95), while the cultivar with the smallest mean value was 'GKB' (22.51).

Identifying Superior Durum Wheat Genotypes Based on the Gluten Index and Minolta b* Values
The applied method of principal component analysis made it possible to investigate the correlations between the examined crop years. The method is used for the analysis of object diversity in regard to quality traits of the winter durum wheat. The most recent method, the GGE biplot model, provided a complete visual evaluation of all aspects of the data by creating a biplot that simultaneously represents both mean performance and stability [31]. The first two principal components, PC1 and PC2, accounted for 93.40% and 3.90% of the GGE sum of squares, respectively, explaining a total of 97.30% of the variation ( Figure 1A) and revealing a differential gluten index performance between the durum wheat genotypes across the years due to the presence of a genotype × environment interaction (G × E). The length of the vectors representing the various years differed, but in all cases, the angle they enclosed was much less than 90 • , showing that the correlation between the years was close and positive. The angles between the vectors and the average environment axis were indicative of whether a year could be regarded as having a 'representative' or 'special' effect. On this basis, 2013 was the clearest example of a 'special' year, while 2007, 2010 and 2018 were located closest to the average. The length of the vectors indicates the extent to which each year made it possible to differentiate between the genotypes. The E14 and E05 vectors were the longest, showing that the durum wheat cultivars could be differentiated to the greatest extent on the basis of their gluten indices in 2005 and 2014. By contrast, the gluten index values of the cultivars differed to the least extent in 2009 (E09). It is clear from the results presented in Table 2, that the highest gluten index values were recorded this year and that the percentage differences between the cultivars were smaller than in the other years. The 'which-won-where' analysis makes it possible to identify the genotype with the most favorable properties for a given environment. The individual sectors represent mega-environments. As all the years were located in the same sector in Figure 1A, and as 'MVP' and 'GKS' were also found in this sector, it can be concluded that all the years had an almost identical effect on the gluten index values of the winter durum wheat cultivars tested; the cultivars with the strongest gluten matrix were 'MVP' and 'GKS', and the cultivar with the weakest gluten matrix was 'MVM', as this was found in the 3rd sector, on the left, at the greatest distance from the environment (year). The analysis of the Minolta b* values again revealed that a single principal component variable was dominant; PC1 and PC2 accounted for 88.10% and 6.20%, respectively, of the GGE sum of squares, explaining a total of 94.30% of the total variance ( Figure 1B  fered, but in all cases, the angle they enclosed was much less than 90°, showing that the correlation between the years was close and positive. The angles between the vectors and the average environment axis were indicative of whether a year could be regarded as having a 'representative' or 'special' effect. On this basis, 2013 was the clearest example of a 'special' year, while 2007, 2010 and 2018 were located closest to the average. The length of the vectors indicates the extent to which each year made it possible to differentiate between the genotypes. The E14 and E05 vectors were the longest, showing that the durum wheat cultivars could be differentiated to the greatest extent on the basis of their gluten indices in 2005 and 2014. By contrast, the gluten index values of the cultivars differed to the least extent in 2009 (E09). It is clear from the results presented in Table 2, that the highest gluten index values were recorded this year and that the percentage differences between the cultivars were smaller than in the other years. The 'which-won-where' analysis makes it possible to identify the genotype with the most favorable properties for a given environment. The individual sectors represent mega-environments. As all the years were located in the same sector in Figure 1A, and as 'MVP' and 'GKS' were also found in this sector, it can be concluded that all the years had an almost identical effect on the gluten index values of the winter durum wheat cultivars tested; the cultivars with the strongest gluten matrix were 'MVP' and 'GKS', and the cultivar with the weakest gluten matrix was 'MVM', as this was found in the 3 rd sector, on the left, at the greatest distance from the environment (year). The analysis of the Minolta b* values again revealed that a single principal component variable was dominant; PC1 and PC2 accounted for 88.10% and 6.20%, respectively, of the GGE sum of squares, explaining a total of 94.30% of the total variance ( Figure  The ranking biplot method was used to evaluate the genotypes on the basis of gluten index and Minolta b* value means and their stability across environments ( Figure 2). The abscissa indicates higher mean values across environments and the ordinate greater variability (poorer stability) in either direction. With regard to gluten index, in all 16 years 'MVP' had the highest mean gluten index, followed by 'GKS', while 'MVM' had the The ranking biplot method was used to evaluate the genotypes on the basis of gluten index and Minolta b* value means and their stability across environments ( Figure 2). The abscissa indicates higher mean values across environments and the ordinate greater variability (poorer stability) in either direction. With regard to gluten index, in all 16 years 'MVP' had the highest mean gluten index, followed by 'GKS', while 'MVM' had the lowest mean value (Figure 2A). Regarding stability, all of the genotypes were almost identical over an average of 16 years. lowest mean value ( Figure 2A). Regarding stability, all of the genotypes were almost identical over an average of 16 years.
Based on the Minolta b* values of the durum wheat cultivars 'MVH' gave the best results, followed by 'MVP', 'MVM', 'GKS' and 'GKB'. 'MVM', which was located almost on the average environment axis (AEA), was extremely stable, averaged over 16 years ( Figure 2B). Interaction between the Years and Winter Durum Wheat Cultivars on the Basis of Gluten Index and Minolta b* Data over 16 Years In the course of selection, breeders need to select simultaneously for several traits, the aim being to develop cultivars with the best combination of these traits. GGE biplot analysis makes it possible to investigate several traits at the same time.
The calculations revealed that two principal component variables had eigenvalues exceeding 1, and these explained 93.3% of the complete variance of the system (Figure 3). The gluten index and Minolta b* vectors formed two separate groups, with an approximately 90° angle between the mean vectors of the groups, indicating the independence of the two traits. The cultivars located at the corners of the rectangle had different combinations of the two technological quality traits tested. 'MVP' was located in the same sector as the gluten index vectors, at one corner of the rectangle, indicating that the gluten index of this cultivar was outstanding. 'MVH' was in a similar relationship with the Minolta b* vectors. 'MVM' was in the sector opposite the gluten index vectors, indicating a weak gluten matrix, while 'GKB' occupied a similarly unfavorable position with respect to the Minolta b* vectors. In the course of selection, breeders need to select simultaneously for several traits, the aim being to develop cultivars with the best combination of these traits. GGE biplot analysis makes it possible to investigate several traits at the same time.
The calculations revealed that two principal component variables had eigenvalues exceeding 1, and these explained 93.3% of the complete variance of the system (Figure 3). The gluten index and Minolta b* vectors formed two separate groups, with an approximately 90 • angle between the mean vectors of the groups, indicating the independence of the two traits. The cultivars located at the corners of the rectangle had different combinations of the two technological quality traits tested. 'MVP' was located in the same sector as the gluten index vectors, at one corner of the rectangle, indicating that the gluten index of this cultivar was outstanding. 'MVH' was in a similar relationship with the Minolta b* vectors. 'MVM' was in the sector opposite the gluten index vectors, indicating a weak gluten matrix, while 'GKB' occupied a similarly unfavorable position with respect to the Minolta b* vectors. Figure 4 allows the durum wheat cultivar with the best approximation to the 'ideal' genotype (in this case, with the best combination of the given traits) to be identified. The value of a genotype is determined by its distance from the theoretical genotype located in the middle of the concentric circles (best for all traits, with perfect stability). The closer the genotype is to this point, the better the combination of desired traits. According to this analysis, 'MVP' came closest to the ideal genotype, with the best combination of these traits. 'GKS' and 'MVH' were at approximately the same distance from the ideal, the former excelling chiefly for gluten strength, the latter for yellow pigment content. 'MVM' and 'GKB' was located at the greatest distance from the ideal genotype.   Figure 4 allows the durum wheat cultivar with the best approximation to the 'ideal' genotype (in this case, with the best combination of the given traits) to be identified. The value of a genotype is determined by its distance from the theoretical genotype located in the middle of the concentric circles (best for all traits, with perfect stability). The closer the genotype is to this point, the better the combination of desired traits. According to this analysis, 'MVP' came closest to the ideal genotype, with the best combination of these traits. 'GKS' and 'MVH' were at approximately the same distance from the ideal, the former excelling chiefly for gluten strength, the latter for yellow pigment content. 'MVM' and 'GKB' was located at the greatest distance from the ideal genotype.

Identification of the Meteorological Factors Determining the Gluten Index and Minolta b* Values
The correlation analysis was performed on a data matrix that provided for 61 meteorological factors (26 rainfall + 23 mean temperature + 12 heat day data) and five cultivars (cultivar mean × 2 traits) ( Table S1). The significant factors for both parameters were

Identification of the Meteorological Factors Determining the Gluten Index and Minolta b* Values
The correlation analysis was performed on a data matrix that provided for 61 meteorological factors (26 rainfall + 23 mean temperature + 12 heat day data) and five cultivars (cultivar mean × 2 traits) ( Table S1). The significant factors for both parameters were shown in Table 3. Table 3. Pearson's correlation coefficients between meteorological factors and the gluten index and Minolta b* values of the five durum wheat genotypes over 16 years. *, ** and *** denote significant differences at the p < 0.05, 0.01 and 0.001 levels of probability.

Meteorological Factor 'GKB' 'GKS' 'MVM' 'MVP' 'MVH' Cultivar Mean
Gluten Index The gluten index data measured for each cultivar and their mean were negatively influenced by the majority of meteorological data with identified effects ( Table 3). The increased mean monthly temperature in June and the mean temperature in the first 10 days of July significantly decreased the gluten index for all the durum wheat cultivars and for the cultivar mean. An increase in the number of heat days in the middle 10 days of May had a positive significant effect on the gluten index of 'MVM' (r = 0.859 ***), 'GKB' (r = 0.533 *), 'MVH' (r = 0.511 *) and on the cultivar mean (r = 0.617 *). The gluten strength of two cultivars ('MVP' and 'MVH') and the cultivar mean significantly decreased in response to greater rainfall in the first 10 days of April, however, 'MVM', which had a weak gluten matrix, responded favorably to the rainfall quantity in the third week of June (r = 0.566 *). Cultivar-specific effects were detected for the mean temperature in the first 10 days of June in 'GKS', and for the quantity of rainfall in November, and in the last 10 days of May and in the whole of May, and the number of heat days in the middle 10 days of June in 'MVP'.
The Minolta b* (MB) value was significantly influenced by seven meteorological factors. The number of heat days in the last 10 days of May intensified the yellow color of all the cultivars and the cultivar mean (r ranged from 0.568 * to 0.777 ***), and a similar effect was observed for the total number of heat days in May, in June and throughout the vegetation period. Rainfall at the beginning of June also had a favorable effect on the intensity of the yellow color (r ranged from 0.546 * to 0.623 **). A cultivar-specific effect was only detected for 'MVH', in which the Minolta b* value was enhanced by the quantity of rainfall in the last 10 days of April, but reduced by the mean temperature in October.

Combined Effect of Meteorological Factors on the Gluten Index and Minolta b* Values of the Cultivars Tested
In the calculations, the rainfall (total rainfall in the months from August to March, quantity of rainfall every 10 days from April to harvest), the mean temperature data (monthly data from October to March, 10-day data from April onwards) and the number of heat days (every 10 days from early May to harvest) were taken as independent variables. The gluten index (gluten strength) was diversely influenced by meteorological factors in the different cultivars. The most accurate regression function could be compiled using only two factors for 'GKS', while four were needed for 'MVM', seven for 'MVP' and 'MVH' and nine for 'GKB'. When calculations were made for the average of the five cultivars, four factors played a decisive role in the development of the gluten index (Table 4).  As several of the winter durum wheat cultivars gave an excellent approximation to the expected value, but the functions were extremely complicated (World S1), data are also given for models where the value of the multiple coefficients of determination (R 2 ) exceeded 0.9. In the latter case the expected value of the gluten index can be estimated to a satisfactory approximation using a simplified function: The Minolta b* value was diversely influenced by meteorological factors in the different cultivars (World S1), except 'MVM' and 'MVP' genotypes. The regression function could be compiled using only one factor for them, as 'No. of heat days in 3rd 10 days of May' played a decisive role in the development of Minolta b* value. Moreover, this factor was determinant in all varieties. The most accurate regression function could be compiled using four factors for 'GKB', five factors for 'GKS' and fifteen for 'MVH'. When calculations were made for the average of the five cultivars, eleven factors played a decisive role in the development of the Minolta b* value (Table 5).  Stepwise regression analysis made it possible to identify meteorological factors with a decisive influence on the genotype-dependent and mean gluten index and Minolta b* value. Averaged over the cultivars, the number of heat days in the second 10 days of May and the mean temperature in the final phase of grain filling were of outstanding importance for the gluten index. Heat days had a positive influence on gluten strength in the early stages of development, and a negative effect in later stages. The negative effect of the precipitation quantity in October can be explained by delays in sowing. Hot weather soon after pollination (No. of heat days in the last 10 days of May) had a positive influence on the Minolta b* value, but in later stages (mean temperature in the first 10 days of June) the effect became negative.
When the paired t-test was used to analyze the similarity of the estimated data (substituted into the gluten index and Minolta b* functions) and the original data set (Table 6), the two sets of data were found to agree at a high level of probability (p > 0.94) for all the durum wheat cultivars and for the cultivar mean. As the five genotypes differ as regards both gluten index and Minolta b* value, it can be concluded from the R 2 values of the regression functions calculated from the cultivar mean for gluten index (0.873) and Minolta b* value (1.000) that the expected values of these important technological quality traits can be reliably estimated on the basis of meteorological factors. The functions for the calculation of the gluten index include data not available until the end of the first 10 days in July, so early estimation is not possible for this trait. The Minolta b* value, on the other hand, is determined by the first ten days in June, so it is possible to predict with considerable certainty 3-4 weeks before harvest.

Discussion
The important objective of the winter durum wheat improvement program is to develop not only heat, drought and frost tolerant, high yielding with better quality but also stable varieties. It was hypothesized that it should be possible to identify meteorological factors with an above-average influence on the gluten index and Minolta b* values of winter durum wheat cultivars. In order to prove this, the correlation between various meteorological factors (monthly total precipitation and mean temperature from August to March, the same factors for each 10-day period from April to harvest, and the number of heat days in each 10-day period in the grain-filling period) and two technological quality parameters was analyzed by means of correlation analysis and stepwise regression analysis.
In the present experiment, the cultivar with the strongest gluten matrix was more stable than those with weak gluten, exhibiting less change in the gluten index in years with warmer weather during grain filling. This, therefore, confirmed in the case of winter durum wheat cultivars the observations made by Ames et al. [16] in spring genotypes. Meteorological factors related to high mean temperature or to the number of heat days are frequent terms in the equations used to predict the mean gluten index of the five cultivars. All of these are preceded by a regression coefficient with a minus sign, i.e., the higher the mean temperature at the end of the grain-filling period and/or the number of heat days during this period, the lower the gluten index and the poorer the gluten strength.
The occurrence of several heat days in the early stages of grain filling (2nd 10 days in May) had a favorable influence on the gluten index, especially in 'MVM', which had a weak gluten matrix, but also in 'GKB' and 'MVH' and the cultivar mean. It can be hypothesized that if heat stress occurs during flowering ('MVM' is a late cultivar, generally heading in the 2nd 10 days of May), it may cause a change in the accumulation dynamics of the various types of storage proteins. Targeted experiments in the phytotron will be required to confirm this. Further studies will also be needed to decide whether the negative correlation between April precipitation quantity and the gluten index of cultivars 'MVP' (−0.733 **) and 'MVH' (−0.600 *) and the mean gluten index of the five cultivars (−0.589 *) was a direct effect or the consequence of other indirect physiological changes.
In all the cultivars, the Minolta b* value was enhanced by the number of heat days in the early stage of grain filling. Furthermore, the values of four of the five durum wheat cultivars and the cultivar mean exhibited a significant correlation with the number of heat days during the vegetation period. These results confirm the findings of Rharrabti et al. [19], who reported that higher temperatures favored an increase in yellow pigment content. The highest Minolta b* value was measured in the present work in samples taken in the driest year (2007: total rainfall quantity in the whole vegetation period: 153.8 mm), suggesting that drought improves the effect of heat stress. Similar results were found in spring durum wheat by Li et al. [21], but it was suggested that intense heat stress might reverse this process, leading to a lower Minolta b* value [22]. This conclusion was confirmed by the observations of Cseuz et al. [25] and Hadži-Tašković Šukalović et al. [26], who did not detect any significant difference between the yellow pigment content in different years under the climatic conditions of the Pannonian Region. In some cases, the year has also been known to influence this trait [24]. In experiments set up between 1991 and 1994 in Szeged (S. Hungary), a substantial difference was recorded between the yellow pigment contents detected in 1993 and 1994, with the highest value in 1994 when the weather was cooler (17.97 • C), with more rainfall (62.6 mm) during the grain-filling period (May and June), with respect to 1993, warmer (mean temperature: 19.38 • C) and slightly drier (56.9 mm). This suggests that cooler weather with more rainfall in the early stage of grain filling had a favorable effect on the yellow pigment content. In Martonvásár (Hungary) durum wheat reaches this phenophase a few days later than in south Hungary, as it is located 200 km north of Szeged. In conclusion, the present data showed the favorable effect of a larger quantity of rainfall in the first 10 days of June, when the early stage of grain filling occurred. The present data thus confirm those of Matuz et al. [24] as to the favorable effect of rainfall during early grain filling.

Conclusions
Analyzing correlations by means of stepwise regression analysis proved to be a satisfactory method for the joint detection of the effect of various factors. GGE biplot analysis on the effect of the year on five winter durum wheat cultivars in 16 growing seasons revealed the years (environments) when selection could best be performed for gluten index and Minolta b*. The cultivars tested differed considerably for both yellow pigment content and gluten strength. The gluten index was outstandingly high for the cultivar 'MVP', followed by 'GKS', which also had excellent stability, whereas 'MVH' was ranked first for Minolta b*, but had poor stability. Although yellow pigment production was determined jointly by a number of meteorological factors, the final value of Minolta b* could be predicted well before full ripening for the majority of durum wheat cultivars. 'MVP' was the best approximation to the ideal genotype (a theoretical cultivar with the optimum combination of the two traits), so on the basis of the long-term database, this cultivar can be recommended for use as a crossing partner in breeding.
As reported by other authors, analyzing correlations by means of stepwise regression analysis proved to be a satisfactory method for the joint detection of the effect of various factors.

Site Description and Set-Up for the Field Experiment
The rain-fed experiment was arranged in a randomized complete block design with three replications and a sowing density of 500 seeds/m 2 . The geo-coordinates of the field were latitude 47 • 30 41.70" N and longitude 18 • 81 83.34" E (Martonvásár, Hungary). The genotypes were sown in the middle of October each year. The soil texture was loam with good P and K supplies and good water permeability. Individual plots measured 0.96 m × 6 m and consisted of eight rows, 12 cm apart. Nutrients were supplied in autumn at a rate of 60:60:60 kg ha −1 N:P:K active ingredients, followed by topdressing with 60 kg ha −1 N in early spring. No fungicide treatment was applied during the vegetation period, and weeds and pests were controlled chemically. The plots were harvested at full maturity using a small-plot combine harvester.
The meteorological data used for the analysis were obtained from the database of an automatic weather station located in Martonvásár, Hungary ( Table 7). The data collected were the precipitation (mm) during the whole vegetation period and specifically during grain filling; mean temperatures ( • C) during the vegetation period and the grainfilling period; the number of heat days (maximum daily temperature ≥30 • C); and the meteorological stress factors characteristic of the vegetation period.

Determination of Gluten Index and Minolta b* Values
Gluten index and b* values of the winter durum wheat samples were determined on semolina produced using a Chopin CD2 mill and purified with a Chopin Laboratory Purifier according to the ICC158 standard (ICC, 1995). The gluten index was determined using a Perten Glutomatic 2200 instrument and a Perten 2015 Centrifuge (Perten Instruments AB, Hägersten, Sweden). The Minolta b* values were determined using first a Minolta CR-300 and from 2016 onwards a Minolta CR-400 chromameter (Minolta Camera Co. Ltd., Osaka, Japan).

Statistical Analysis
Analysis of variance was performed for years and cultivars using the one-factor random block design model combined with the year in the MSTAT-C software package (Michigan State University, East Lansing, MI, USA). The significance of differences between the mean values of the main factors was determined using Tukey's HSD test. Complex correlations between the combined effects of several factors were analyzed by means of stepwise regression analysis. Interactions and correlations between years and cultivars were determined using principal component analysis, while the best-performing genotype was selected by means of principal component biplot analysis (GGEbiplot Pattern Explorer, Version 8.1) [31].
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/plants11010113/s1, Table S1. Pearson's correlation coefficients between meteorological factors and the gluten index and Minolta b* values of the five durum wheat genotypes over 16 years.